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We consider a model of inflation consisting a single fluid with a time-dependent equation of state. 
In this phenomenological picture, two periods of inflation are separated by an intermediate non- 
inflationary stage which can be either a radiation dominated, matter dominated or kinetic energy 
dominated universe, respectively, with the equation of state w = 1/3, or 1. We consider the toy 
model in which the change in w happens instantaneously. Depending on whether the mode of interest 
, leaves the horizon before or after or between the phase transitions, the curvature power spectrum 

can have non-trivial sinusoidal modulations. This can have interesting observational implications 
(~^) for CMB anisotropies and for primordial black-hole formation. 

(N' 

I. INTRODUCTION 

Inflation has emerged as the leading paradigm for early universe cosmology and structure formation. Basics predic- 
tions of inflation are in good agreement with cosmological observations. Namely, simplest models of inflation predict 
almost scale invariant, almost Gaussian and almost adiabatic fluctuations on cosmic microwave background (CMB) 
which are accurately measured in recent cosmological observations [l[. Nonetheless, it is interesting to consider more 
elaborate models of inflation which can predict observable deviations from these simple predictions. In particular, 
models of inflation with local features may be interesting. Observationally, these models are employed to address 
the glitches in the CMB angular power spectrum on scales £ ~ 20 — 40. Theoretically, one can construct different 
scenarios which can contain local features Models with local features may originate from high energy physics, 

£C) . particle creations, field annihilations, change in sound speed or time variations of the Newton constant during in- 
^sO 1 nation [17H3ll |. Many of these models are based on multi-field or multi-fluid scenarios. As a consequence, there are 
C*") , always iso-curvature perturbations which may be constrained from CMB observations. 

■ In this work we consider a phenomenological model with multiple inflationary stages. Different stages of inflation are 
separated by an intermediate non-inflationary period. In our model, these multiple inflationary stages are realized by 
changes in the equation of state w for a single fluid. During inflation w ~ — 1 while in the intermediate non-inflationary 
stage we have 1 + Zw > 0. Particular interests are the cases in which the intermediate non-inflationary stage has the 
equation of state u> = u>2 = l/3,0orl, corresponding respectively to a radiation, matter or kinetic energy dominated 
universe. Having this said, we should emphasis that this is a phenomenological study and a dynamical mechanism 
causing the jump in w has yet to be constructed. We shall briefly present a simple scalar field model which can 
provide a simple dynamical mechanism for changing w. Idea similar to this line of thought was studied in 32] in the 
context of MSSM inflation 

The rest of the paper is organized as follows. In section |TT] we present our setup and background equations. 
In section IIIII we present the general perturbation equations with appropriate matching conditions. The resulting 
transfer function of the outgoing perturbations for arbitrary u>2 7^ is given in Section HVl The special case of W2 = 
is considered in section [V] The conclusion and discussions are given in section [VT] and some technical issues are 
relegated to appendices. 
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II. THE BACKGROUND 



Here we present the background evolution of a universe filled with a single perfect fluid with an arbitrary but 
constant w — P/ p where p and P are the energy density and pressure, respectively. The background space-time is 
assumed to be a flat FLRW universe, 

ds 2 = -dt 2 + a{t) 2 dx 2 = a 2 (j])(-dr) 2 + dx 2 ) , (1) 

where 77 defined by dr\ = dt/a(t) is the conformal time. 

To be specific we have the following picture in mind. We have three distinct stages of an expanding background in 
which two inflationary periods are separated by an intermediate non-inflationary stage. The first inflationary stage 
continues till 77 = 7712 and during this period the fluid driving inflation has a constant equation of state w = w±. In 
order to support inflation we require 1 + 3wi < 0. We assume that at 77 = 7712 the equation of state changes sharply 
from w\ to W2 such that 1 + 3w2 > and the first stage of inflation is terminated. As specific examples we shall 
consider the important cases of W2 — 1/3, and 1, corresponding respectively to radiation, matter and kinetic energy 
dominated universes. The third expanding stage starts at 77 = 7723 when w goes a second abrupt change from W2 to 
W3. In order to support the final stage of inflation we assume that 1 + Sw^ < 0. The time when the inflation ends is 
set to be 77 = rj e = followed by a (p)reheating era. In summary, we have two inflationary stages with equations of 
state w = wi and W3 separated by an intermediate non- inflationary stage with w = w%. 

One may wonder how dynamically these jumps in w can be realized in a consistent way. In Discussions Section we 
present a simple scalar field model which can mime this behavior. However, in this and the following sections, we 
shall proceed phenomenologically assuming that there exists a dynamical mechanism which can cause these changes 
in w. For our analytical analysis we proceed with the arbitrary sharp changes in w. We note that physically it is 
expected that the process in which w undergoes large changes will take some finite lapse of time. Therefore we also 
consider numerically the case when there is a short but finite duration of the phase transitions. We shall also compare 
our analytical results with sudden changes in w to those obtained numerically in which the change in w takes a finite 
lapse of time. 

With this picture in mind, now we present the background equations. Using the energy conservation equation and 
denoting the initial conditions with subscript 0, the evolution of energy density is given by 

( v -3(1-H») 

P = (*>(%) ■ ( 2 ) 

Here and below the subscript collectively denotes the time of phase transitions, so it corresponds to either 7712 or 
7723 depending on which period is studied (see Eq. ([7]) below for further details). The Friedmann equation for a flat 
universe is 

3Mpi 2 W 2 = a 2 p, (3) 

where % = a' /a is the conformal Hubble parameter and the prime denotes the derivative with respect to conformal 
time 77. The Friedmann equation for w 7^ — — can be integrated to 



where ao = 0(770) and 



a(r7)=«o(y(r7-77o) + l) , (4) 



(5) 



Taking the conformal time derivative of Eq. ((4]), one obtains 



^ = —7^ • (6) 



It is easy to see that both the scale factor 0(77) and the conformal Hubble parameter H must be continuous at the 
time of phase transition 77 — 770 when w undergoes a sudden change. 
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We label quantities at the three stages with 1, 2 and 3. For example "Hi (77) is the conformal Hubble parameter 
during the first stage whereas H12 is the value of at the time of first phase transition 77 — 77121 H12 — ^( 7 7i2)i 
and so on. As a result one has 

( u„„ 

for 77 < 7712 , 

for 7712 < 77 < 7723 , (7) 
for 77 > 7723 , 



nv) = { 





1- 


h?M77- 






H12 




1 - 


h Hliiv- 






H23 




1 - 




f?23)/% 



in which /3, are defined as in Eq. ([5]) with w replaced by u>i of each stage. Here we assume that the first stage of 
inflation starts at 77 — > —00 and the second stage of inflation ends at 77 e = 0. Note that the Hubble parameters at two 
phase transitions are related by 



7^23 — 



12 



1 + ^12(7723 - Vl2)/&2 

while the Hubble parameter at the end of inflation is given by 

%3 



1 - H23 7723/^3 



(8) 



(9) 



Note that, as long as the intermediate non- inflationary stage corresponds to a universe dominated by an ordinary 
matter (w > 0) one has 



fti2 > rl 



23 



23 ■ 



(10) 



Alternatively, it may be useful to work with the number of e-folds as the clock dn = Hdt = T-Ldrj. With the scale 
factor given by Eq. (j4]) one obtains 



71(77) =n + p In (1 + p-^ri - Vo)) ■ 



(11) 



In particular, the number of e-folds of the second non-inflationary stage AA^ = 77.23 — ^12, where 7123 = 77,(7/23) and 
7ii2 = 77(7712), and that of the third inflationary stage AA3 = n e — 7723, where n e = n(j] e ) = n(0), we obtain 



57^ = 1 + ^2 ^lafej - V12) = < 
ri-23 



tie 



= e AN 3 //3 3 



(12) 



III. THE PERTURBATIONS 



In this section we study the perturbation equations in details. We study the behaviors of the comoving curvature 
perturbations TZ or the Bardeen potential $ which are gauge invariant. Some technical details are described in 
Appendix [A] For a review see, e.g. (33|. 

For a universe filled with a single fluid with the known equation of state parameter w and sound speed c s one 
obtains the following equation for the Fourier space mode function of the comoving curvature perturbation: 

Tl'i + + C 2fc 2 ^ k = , (13) 

where 

z = a(7?)M P1 V3(l + w)/c s . (14) 

Note that c s is defined as 8P C — <? s 5p c in which the subscript c indicates that the corresponding quantities are measured 
on the comoving hypersurface (on which the fluid 4- velocity coincides with the unit normal to the hypersurface). 

One can easily solve Eq. (|13l) in each phase with constant values of w and c s . At the time of transition, we need 
two matching conditions in order to match the outgoing solutions to the incoming solutions. 
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The first matching condition is the continuity of the curvature perturbation itself, 

[ftk]t =0, 



(15) 



where [X]~t denotes the difference in the value of quantity X after and before the transition: = X(r/ + ) — X(r)-). 

Geometrically, the continuity of 7Z can be interpreted as the continuity of the extrinsic and intrinsic curvature at the 
three-dimensional spatial hyper-surfaces located at r/ = 7712 and r/ = 7723, to be consistent with the Bianchi indentity. 
We also need another matching condition for the time derivative of 1Z. Note that the Eq. (jT3l can be rewritten by 



(16) 



— ( ^-(1 + w)TZ' k ) + a 2 (l + w)k 2 1Z k = . 



By integrating the above equation in a small range around the phase transition, the last term vanishes and one obtains 
the second matching condition by 



1 + M V 



= 0. 



(17) 



Alternatively, one can obtain the matching condition (|17|) in a different way. From the continuity of the extrinsic 
and intrinsic three-dimensional hyper-surface at the time of phase transition, we also conclude that the curvature 
perturbation on the shear-free hypersurface (Newton gauge) $ is continuous across the transition surface, 



[$]+ =0. 

Then by noting the relations between $ and 1Z (see Appendix [A] for details ) 

2$' + (5 + 3w)H<i> 



K 



and 



m{i + w) 



2c?fc 2 



(18) 



(19) 



(20) 



we see from Eq. ([2H)l that the matching condition (|T7|) implies the continuity of $. 

Now we solve the equation of motion for TZ^. For a constant w and c s , we have z 1 / z = a' /a = %. Hence Eq. (fT3 
simplifies to 



dx 2 



x ax 



Ku = 0. 



where 



x = — c s k{rj - 7] + (3H x ) . 



(21) 



(22) 



The solution is given by 



CiHW (x) + D.H^ (x) 



3(w-l) 
2(3w + 1) 



(23) 



where C\ and D\ are constant of integrations, and Hu (x) and Hu (%) are the Hankel functions of the first and 
second kinds, respectively. Note that during inflation f3 < (for slow-roll inflation (3 ~ —1) and the above general 
definition of x yields 



x = xi(rj) = -c sl k(r) - i] 12 + PiH^) , 



(24) 



during the first period of inflation. 

Our goal is to find the curvature perturbation 7Z k at the end of inflation 77 = 0. The power spectrum V-jz is defined 

by 



(KtKv) = (2nfP n (k) S 3 (k + k') , Vn = ; Pn{k) = |^ k | 2 



(25) 
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where IZ^ is the normalized positive frequency mode function. At sufficiently early times, the solution should approach 
the Minkowski positive frequency mode function. That is, for r\ — > — oo, 

g— ic s kr\ 

Imposing this initial condition on the solution (1231) and using the asymptotic form of the Hankel function given by 
Eq. (|B1[) . we find at the first stage of inflation, D\ = and 

n^ri) = C 1 x l {r,) v -H^{x l {r l ))-, V < m , (27) 

where 

2Mpia(r]i2) V3fc(l + wi)/ 

Considering the slow- roll limit in which wi = — 1 + 2ei/3 and z/i ~ 3/2 the above equation results in the following 
power spectrum at r\ = rji2 for the modes which leave the horizon during the first stage of inflation: 



Hi 



PizM^-y—^ , (29) 

8tHA/ P i c sl e 1 

where the slow-roll parameter e is defined by 

e=~ = \{l + w). (30) 

For those modes that remain superhorizon until the end of inflation, the curvature perturbation is conserved, and we 
have 

Pn(Ve) = Vizivu) * 2 ,f k • (31) 

This is the standard results for single field inflation [3~ij ]. 

Applying the general solution (f23| to the second and third stages, we have 

ftkfa) = C< 2 x?H£> (x 2 (ri)) + D 2 xpH™ (x 2 ( V )) ; Vl2 < V < V23 , (32) 
ftkfa) = C 3 a%>Hg> {x 3 (r,)) + D 3 x?H%> (x 3 ( V )) ; ry 23 < r, , (33) 
where Xi(r]) (i — 2, 3) are defined in accordance with the general definition ([22]) . 



x 2 = c s2 k(r) - r)i2 + foHw) ; v 2 = \- fi 2 = 3 / W2 ^ 

2 2(3w 2 + l) 



a; 3 = c s3 fc(r, - V23 + p 3 U^) ; v 3 = - - /3 3 = , . (34) 



1 _ 3(^3-1) 

2 P3 2(3w 3 + 1) 

Here we assume that w 2 = c~ 2 7^ 0, so the intermediate non-inflationary stage is not a matter-dominated universe. 
The case when the intermediate stage is matter dominated with w = c? s2 — is considered separately in Section |Vl 

As usual, we are interested in modes which are super- horizon at the end of inflation x 3 {rj e ) <C 1 where r\ e — > 0. 
Using the asymptotic form of the Hankel function, we obtain 

i 2 V3 

ft k (r?^0)~ T{v 3 ){C 3 -D 3 ). (35) 

7T 

It is useful to define the transfer function T for the power spectrum as 

Vn(v = 0) = TV n ,(v = 0), (36) 

where Viz-i {f} — 0) is the power spectrum at the end of inflation if there were no transition and w = w\ throughout 
the inflationary stage, as calculated in Eq. ([3"Tj). Assuming v\ ~ v 3 ~ 3/2 we obtain 

1 - Idp • {61) 

Thus any non-trivial effect due to change in w is captured by a non-trivial transfer function T ^ 1. The details of the 
calculation of the coefficients C 2 , D 2 , C 3 and D 3 are given in Appendix [Cj 
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FIG. 1: Here we plot T-i/c a (left) and the dimensionless number H/kc s (right) as a function of the number of e-folds n. As 
usual, during inflation the ratio H/c s increases while during the second stage (non- inflationary period) it decreases. In the left 
plot, the upper, middle and lower horizontal lines correspond, respectively, to k = 5fci2, k = 0.5fci2 and k = 0.05fci2. As can be 
seen in both plots these modes enter the horizon at different times. 



IV. TRANSFER FUNCTION 



In this section we calculate the transfer function T which encodes the effects of change inw. As mentioned above, in 
this section we assume that the intermediate non-inflationary stage is not matter dominated so W2 = c 2 s2 7^ and our 
formulas (ICip - (|C4l) are valid. The case in which W2 = c 2 s2 = is studied separately in Section [V] Also we provide the 
general formula for arbitrary w 2 7^ and then consider the particular example 102 = 1/3 and w 2 = 1 corresponding, 
respectively, to a radiation dominated and kinetic energy dominated universe. 

The behaviors of the power spectrum depends on the ratio kT-L/c s . Depending on this value different situations 
arise. To be specific let us define 

ki2 = -5 , fc 23 = . (38) 

PlCsl P3C S 3 

Here k\ 2 represents the mode which leaves the horizon at the time of the first phase transition whereas &23 is the 
mode which leaves the horizon at time of the second phase transition. Note that k\ 2 > ^23 as long as the second stage 
is non-inflationary. With the above definitions of the characteristic wave numbers, there are three categories of the 
modes. The behavior of H/c s as well as that of the ratio of T-L/c s to k for three typical values of k is depicted as a 
function of the number of e-folds in Fig. [TJ 

The first category contains the modes k > fci 2 , which remain sub- horizon until the second stage of inflation, For 
this category we have £1(7712), £2(^12), £2(7723), £3(7723) > 1. 

The second category contains the modes ^23 < k < ki2, which leave the horizon during the first inflationary 
stage, re-enter the horizon during the second non-inflationary stage, and finally exit the horizon during the second 
inflationary stage. For this category we have £1(7712), £2(7712) "C 1 and £2(7723), £3(7723) 3> 1. 

The third category contain the modes k < k 2 ^, which leave the horizon during the first stage and always remain 
super-horizon until the end of inflation. For this category we have £1(7712), £2(^12), £2(^23), £3(^23) *C 1. 

Now we study each category of the modes in turn. Let us start with the first category, modes which remain 
subhorizon until the second phase of inflation, k > k\ 2 (for which £1(7712), £2(7712), £2(^23), £3(^23) 1)- The 
coefficients C2, D 2 , C3 and D3 are calculated in Appendix O The amplitude of the curvature perturbation at the end 
of inflation is obtained as 

17^(77 = 0)1 - /23^3-h>3) *f 12 ir: 

n £3(?723) 3 ~ 2 \X2{Vl2), 

x Vl + sin (2£ 3 (1723) - 7^3) sin( x 2 (7712) - £2 (7723)) • (39) 
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Then the transfer function T is calculated to be 

I*/* ~ V2f 23 ^^^ f^ip^Y 2 2 ^ + sin (2a;3(rta) _ ^ sin^r^) - £2(7723)) 

fc V 1 " 4 / * \ 4-1/3 (Pxcsikx^^- 1 



^2/23 h— h— 



&12/ \ fc 23/ \/?3C s 3fe3 

/3 2 c S 2 k 
hc s i ki2 




•kvz sin 



/33C S 3 &23 



(40) 



where we have defined 



sgn(l + 3wj) (1 + Wj) c S j 
sgn(l + 3i0j) (f + Wj) c s 



in which sgn(x) is the sign function; sgn(x) = +1 (—1) for x > (x < 0). Note that in the continuous limit where 
wt — > Wj we obtain the expected result that — > 1 corresponding to no sharp transition. Since a change in w 
naturally causes a change in c s too, it is the combination /y- which controls whether or not we have a non-trivial 
phase transition. 

Now consider the second category of the modes &23 < k < ki2 , which leave the horizon at the first stage of inflation, 
re-enter during the intermediate stage and cross the horizon again during the second stage of inflation (for which 
£1(7712), £2(7712) <C 1 while £2(7723), £3(7723) 3> 1). From the result given in Appendix [Cl we obtain 

\K k (0)\ ~ ^/232^- 1/2+ ^r(^ 1 )r(-^ + f) S 3(7723)-^ + ^£2 (7723) l/2 -5 



x y/(l - shi(2x 2 (?723) - 3^)) (f + sin^fe) - ™ 3 )) • (42) 



The transfer function is given by 



T 1/2 ~ /23r(l-^ 2 )£3(7723)~' y3+2 £2(7723) ,y2 " 5 

7T 

x\/(l - sin (2^2(7723) - 3^ 2 tt)) (f + sin (2x3(7723) - nv 3 )) 

v*-" 2 , rn J k y V3+ i fp 2Cs2 k 
^ /23r(i - u 2 ) 7— 

7T \K23 



-2(3 2 C S 2 k 

fac s3 k 



2:1 




3i/ 2 tt ] 4/I + sin I v 3 -k . (43) 



Finally, consider the third category, k < k 23 , corresponding to the modes which leave the horizon during the first 
stage of inflation and never re-enter the horizon until the end of inflation. Since the curvature perturbation is conserved 
on superhorizon scales, the amplitude of these modes are simply given by the standard result given in Eq. (|31[) . and 
the transfer function is trivial; T ~ 1 . 

The results obtained above show that neither the amplitude nor the scale-dependence of the power spectrum is 
the same as the standard case for the modes k > k 23 - A typical example of the transfer function is shown in Fig. [5] 
The power spectrum is highly oscillatory as a function of momentum. This non-trivial behavior is a result of the 
scattering of the initial wave function by the two phase transitions which lead to a mixing of the negative frequency 
modes, which are absent initially. As a result the subhorizon modes at the second phase of inflation are no longer 
purely positive frequency. 

We note that in the limit k 3> &i2 we obtain the scaling property T oc k 2vi ~ 2u3 from Eq. (|4U|) . The non-decaying 
sinusoidal modulation on top of this mild scale-dependence is because of the assumption that the changes in w take 
place abruptly. Under this assumption the small scale modes, no matter how deep inside the horizon they are, are all 
affected. However, if we allow a finite time-scale for the change in w, say At — S, then the sinusoidal modulations on 
the power spectrum die out for frequencies bigger than 5 _1 and the power spectrum reaches its almost scale-invariant 
value at the end of inflation. This behavior is seen in Fig. [3j which is obtained numerically for an example of smooth 
changes in w. 

Before closing this section, let us explore the dependence of the power spectrum enhancement as a function of 
the duration of the intermediate non-inflationary stage, AN 2 - Here we propose two methods to see this behavior in 
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FIG. 2: The transfer function T is plotted for the case where 102 = 1/3, i.e. the non- inflationary stage is radiation dominated. 
The numerical parameters are wi = —0.99, 103 = —0.93, c s i = c S 3 = 1, c S 2 = v^2- Three distinct categories are recognized 
as discussed below Eq. (|38p . corresponding to k > fci2, ^23 < k < ki2 and fc < fc23. The blue curve is from our analytical 
solution (|37|l while the red dots are from the full numerical result. 




FIG. 3: The transfer function T (the red dashed line) for the same situation as in Fig. [2] above but with a relatively mild phase 
transition in which the change in w takes place in 1/3 of an e-fold. As expected, for sufficiently small scales the sinusoidal 
modulations disappear. The blue solid curve is the analytic expression for the sharp phase transition. 
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FIG. 4: The transfer function T is plotted for the case where u>2 = 1, i.e. the intermediate non-inflationary stage is a kinetic 
energy dominated universe. The other numerical parameters are the same as in Fig. [2] with the same physical interpretations 
for the modes behaviors. 

which each has its own advantages. We especially concentrate on Eq. (|43[) which to good approximation shows the 
behavior we are looking for. Firstly note that for a fixed value of T-L e and assuming ~ — 1, the value of fc 2 3 is nearly 
determine by AiV 2 and as a result we would like to compare the values of T for different AiV 2 at the mode k — fc 23 . 
This is the extreme limit of the validity of Eq. (|43p which holds for the modes in the range fc 2 3 <C k <C fci 2 . Obviously, 
one should have fc 2 3 <C fci 2 if one uses the above approximation, which is the case when A7V 2 is sufficiently large. 

The first approach is to set the parameters w\ and W3 to constant values and look for the AiV 2 dependence in the 
transfer function. An inspection of Eq. (j43|) shows that for the scale k — fc 2 3, there is no AiV 2 dependence on the 
amplitude of the transfer function. As a result, the transfer function will not vary much as a function of AiV 2 . This 
argument is supported by Fig. [5l in which the value of the transfer function at k = fc 2 3 is shown as a function of AA^ 2 . 

The second approach is to set the power spectrum of the two inflationary stages equal to each other and vary A./V 2 . 
In this case we impose the condition, 

U ^ ~ ^23 . (44) 

eiC s l £3C S 3 

As a result, we should change e.g. £3 for fixed values of ei and c s i, when we vary AN2, since the Hubble parameter 
changes considerably during the intermediate stage. In fact, using the above equality, one has 

e 3 ~ £ie -2AJV 2 (i+i//3 2 ) (45) 

Noting that 1 + 103 = 2e3/3, and / 2 3 oc 1/(1 + W3), one can conclude from Eq. P3")l that the transfer function behaves 

the enhancement is exponentially larger for larger values of A7V 2 , which is a very interesting 
phenomenon. This feature is also supported from the exact numerical result shown in Fig. [51 

V. MATTER-DOMINATED INTERMEDIATE STAGE 

The analysis in the previous sections are valid as long as the sound speed and equation of state in each stage do not 
vanish. However, if the intermediate stage is a matter-dominated universe then the previous results are not applicable 
since the matching conditions as well as the equation of motion of curvature perturbation are singular. In this section 
we obtain the power spectrum of curvature perturbation for this case. 



10 




11 



In this special case, it is better to work with the curvature perturbation on the Newtonian slice $ (i.e., the so-called 
Bardeen potential), since the equations behave properly when written in terms of $. The relations between 1Z and $ 
are given in Eqs. (fl~9f and (|20[) . Eliminating 1Z from these formulas results in the following second order differential 
equation for 



*" + 3K(1 + 4)$' + [ c 2 s k 2 - 3(w - c 2 jH 2 ]<S> = , 



(46) 



where 



EL 
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3{l + w)H 



(47) 



It is crucial to realize that c w is not the same as c s in general. Note that c s is defined as SP C = c 2 s 8p c where SP C and 
8p c are the adiabatic pressure and energy density perturbations on comoving slices. Only for a universe dominated 
by a perfect fluid, we have c w — c s . 

Note that the equation of motion for the Newtonian potential is well-defined even for the case in which both w 
and c s vanish. In order to solve the above equation we need two matching conditions for $ and As before, the 
continuity of the intrinsic and extrinsic curvatures at the surface of phase transition implies that $ is continuous 
as given in Eq. (TT51) . On the other hand, using Eq. ([TO)) and noting that both TZ and $ are continuous across the 
transition one obtains the second matching condition, 



(1 + w)H (1 + w) 



0. 



(48) 



Now we have the necessary information to obtain the solution for $. At the first stage, the solution is 

* 1 =A l x?- 1 H%l 1 (x 1 ). 
with the same definition of xi(r]) and V\ as it is in Eq. (|22|) . Using the relations (l20l) and (|27|) one has 



(49) 



--13,(1 + w 1 )C 1 . 



(50) 



For the second stage with w-i = c S 2 = , one has 



$9 = Ao 



H12 



(r) - V12) + 1 



B 2 



(51) 



Finally, for the third stage one has 



$3 = x 



AzHllUx^ + BzHllUxs) 



(52) 



Since we are interested in the power spectrum of the comoving curvature perturbation at the final stage, we use 
Eq. ([Tl]) to obtain 



n k = — 



2 x 3 3 



3(1 + w 3 )/3 3 
From this the transfer function is obtained as 



A 3 Hil\x 3 ) + B 3 H^(x 3 ) 



(53) 



rpl/2 = 



21 A, - B« 



31^1(1 + ^3)1^ 



(54) 



The explicit expressions of A3 and B3 in terms of Ci are computed in Appendix [Cl 

In order to obtain an approximate expression for the transfer function, firstly note that (%23/%i2) 5 *C 1 for u>2 > 0. 
Besides that all of Pi and Wi are of the order of unity. Using this information as well as the asymptotic behavior of 
the Hankel functions, we find for the first category, fc > ki2, 
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FIG. 7: The same plot as in [2] but with W2 ~ 0. 



and for the second category, k 2 3 < k < k\ 2 
1/2 ^ 2^- 2 ^(5 + 3^) 



50? #3(1 + 103) 



r>i - 1) 



k 

k 23 



-1/3+3/2 



1 + sin 



2k 

k 23 



&23 < k < ki 



(56) 



In Fig. [7] both numerical and analytic results are shown, which are in very good agreement with each other. 
Interestingly, the power spectrum on small scales is highly enhanced due to the sharp phase transition to the matter- 
dominated era. This agrees with the result obtained in [35| , and will have interesting implications for the primordial 
black hole formation 136 391. 



VI. DISCUSSIONS AND CONCLUSIONS 



In this work we considered a universe in which inflation has multiple stages, separated by intermediate non- 
inflationary stages. To be specific, we considered two inflationary stages separated by an intermediate non-inflationary 
universe. To simplify the calculation, we also assumed sharp transitions between the stages, and studied the case 
when the intermediate stage is either radiation dominated, matter dominated or kinetic energy dominated. 

As well known, in the absence of isocurvature or entropy perturbations, the comoving curvature perturbations are 
frozen on supcrhorizon scales, TZ = 0. So the transitions do not affect supcrhorizon modes. However, sub-horizon 
modes are significantly affected. Hence the final amplitude of the comoving curvature perturbation depends crucially 
on the equation of state parameter of the intermediate stage, w — w 2 , as well as on the time when the mode of interest 
leaves the horizon. 

For modes which leave the horizon at the second stage of inflation, we see a sinusoidal modulations. If the phase 
transition is arbitrarily sharp, then these sinusoidal modulations persist down to infinitely small scales. This is an 
artifact of our assumption that the change in w is instantaneous. In a realistic situation in which the change in w 
takes some finite time scale, then on sufficiently small scales the power spectrum retains it pure vacuum form with 
no sinusoidal modulations. 

We found that the maximum enhancement in the power spectrum occurs for modes which leave the horizon either 
at the time of the first phase transition k — k 23 or at the time of the second phase transition k = ki2- For the case 
of a matter-dominated intermediate stage, w 2 = 0, we found that the power spectrum increases towards small scales. 
This can have interesting implications for the primordial black-hole formation. It would be very interesting to study 
the possibility of the over-production of primordial black-holes in our model. 
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FIG. 8: The numerical plots for equation of state (Left) and the behavior of scalar field (Right) for the model in (|57[) . 

In our analysis so far we have not provided a dynamical mechanism which causes a rapid change in the equations 
of state. Here it may be worth mentioning a simple toy model which can mimic this situation. Consider a single 
scalar field with the canonical kinetic term which is slowly rolling down the potential. This corresponds to the first 
inflationary stage. When the inflaton field reaches near the minimum of the potential and starts to oscillate, the 
intermediate non-inflationary stage commences. The equation of state of this stage depends on the shape of the 
potential near its minimum. For a simple quadratic potential, the universe behaves as matter dominated. 

Now we add a small cosmological constant to the potential which would not affect the dynamics of the first 
inflationary stage. After the energy density of the scalar field is diluted sufficiently due to damped oscillations, the 
cosmological constant starts to dominate and the second stage of inflation begins. 

Note that this is a single field scenario, so our previous analysis are applicable here. However in order to terminate 
inflation, one may need an auxiliary field. Nevertheless, one can assume that this field is sufficiently heavy so that it 
does not affect the large scale (CMB scale) perturbations, as in the hybrid inflation scenario [40ll4lT|. 

The above model can be realized by the potential as simple as 

V = ^m 2 4> 2 + V Q . (57) 

We have used this for the numerical plots in Figs. |S] and O As it is clear from the figure, the first phase transition can 
be quite sharp, while the second transition turns out to be relatively mild. Note that this case is different from what 
we have plotted in Fig. not only because of the mild second phase transition but also because the sound speed is 
equal to unity throughout all the stages, even for the intermediate stage. This is because for the scalar field with the 
canonical kinetic term the sound speed is equal to unity. 

If we relax the assumptions we adopted in the present paper, and consider a multi-field or multi-fluid system, then 
we may be able to construct more realistic models. However the price to pay is the complication of the equations of 
motion and the possible appearance of significant entropy /isocurvature perturbations which are severely constrained 
by the current observational data. Nevertheless it would be interesting to look for such models. 
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FIG. 9: Two numeric result for two different phase transitions. The solid blue line corresponds to a sharp phase transition 
while for the dashed red line the second transition is mild. In both cases we set c s i = c a 2 = c S 3 = 1. The other parameters are 
the same as in ((2| 



Appendix A: Relations between $ and 1Z 



In this appendix we obtain the relations between $ and 1Z as given in Eqs. (JTHJ) and (|20l) . The scalar perturbation 
of the metric is 

ds 2 = -(1 + 2A)dt 2 + 2ad l Bdx l dt + a 2 [(l + 2F)5 ij + 2d l d j E\dx l dx\ (Al) 
The 00, ii, Oi and i =/= j components of Einstein equations are 

3H(F-HA) + ^[F-m] =^ (A2) 

SP 

(2H + 3H 2 )A + HA — F — 3HF = T (A3) 

v ; 2M 2 p y ' 



—F + HA =-P^- v (A4) 

2M 2 V ' 

8 + H6 - A- F =0 (A5) 

in which 8p and 5P are the energy density and the pressure perturbations. We have defined 

9 = a 2 (E-B/a), (A6) 

and v is the fluid velocity scalar potential defined as u,: = div in which is the fluid four-velocity vector. A 
hypersurface on which v — is called the comoving slice. 

The gauge invariant variables Newtonian potential ^ and the Bardeen Potential $ are given by 

* = A - 6 , $ = F - HO. (A7) 
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Furthermore, the curvature perturbation on constant energy hypersurface £ and the curvature perturbation on co- 
moving hypersurface 1Z are given by 

TT 

C = F Sp TZ = F + Hv. (A8) 

P 



We work in an isotropic background with no anisotropic pressure perturbation, so from Eq. (|A5I) one concludes that 

To find Eqs. lfT9j) and (|20|) it is very helpful to go to the comoving slice on which v = and F = 7Z. In this gauge 
from Eq. (|A4I) we obtain A = 1Z/H. Plugging this into Eq. (| A7|) to remove 9 and using $ = — \P yields 

TT , , 

K = $ - (i + H<Pj . (A9) 

Using H = —3(1 + w)H 2 jl and going to the conformal time one can easily check that Eq. (IA9|) is equivalent to 
Eq. (fT9|) . To obtain Eq. (|20)) we have to use the remaining two Einstein equations, Eqs. (|A2I) and (|A3|) . The key 
point is that the sound speed of perturbations c s is defined on the comoving slice as 5P C = c 2 s 8p c where SP C and Sp c 
are the pressure and energy density perturbations on the comoving slice. Using this relation in Eqs. (|A2I) and (|A3[) 
and noting that A = TZ/H on this slice, one obtains 

*=|Sf R - (AID) 

As above, using H = —3(1 + w)H 2 /2 and going to the conformal time one can easily check the equivalence between 
Eq. (fATOj) and Eq. (|2"0)). 

Appendix B: Asymptotic behavior of the Hankel function 

Here we recapitulate the asymptotic behavior of the Hankel functions. For large values of the argument one has 



HP (Z » 1) ~ J J- e«( z — ' / 2 ) , (Bl) 
" ittZ 



whereas for small values of the argument, 



H^(Z<1) ~ -r(H)(^J (-) e-*f("-IH). (B2) 



Also the following identities are helpful: 



and 



JI£(Z) (B3) 



H?XZ)HI*>\Z) H?>\Z)H?\Z) = -4 ■ ( B4 ) 

Appendix C: Transfer Function for Outgoing Solutions 

By solving the matching conditions here we calculate the coefficients of outgoing solutions C2,D 2 ,C3 and D 3 for 
u)2 7^ and A2,B 2 ,A 3 and B 3 for W2—O in terms of C\. We consider each case separately. 

1. TO2 / 



Consider the case where the intermediate non-inflationary stage is not matter dominated, w 2 — c? s2 7^ 0. The 
general solution for the comoving curvature perturbation at the intermediate non-inflationary and second inflationary 
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stages are given in Eqs. (l32j) and (l33l) . respectively. With the help of the identities (|B3I) and (|B4[) and imposing the 
matching conditions (| 15[) and (|17p at 7712 and 7723, we obtain 



C 2 - 
D 2 = 

C 3 = 



4i x 2 (?7i2) I/2 



4^2(7712)^ 
tt £2(7723)'' 



-C\ H^{x 1 {r ]l2 ))H^_ 1 {x 2 { m2 )) hMlUx^^H^ix^)) 



C 2 (HllHx 2 ( m3 ))Hi*Ux 3 ( m3 )) - hsHlll^i^Hifix^)) 



Do = 



4z xs^sY 3 - 1 

+ D 2 (H^(x 2 (r l23 ))Hlf_ l (x 3 ( m3 )) - f 2 sH^ ) _ 1 (x 2 (r l23 ))H^(x 3 ( m3 ) 

^£2(7723)" 



4i x 3 (7723) ly3_:L 



C 2 (^W^^lMto)) - /23< 1) _ 1 (X 2 (?723))^ ) (^3(^3)) 

+ D 2 (H^^im^H^ix^)) - /23^ ) _ 1 (x 2 (r/23)X 1) (^(r/23) 



(CI) 
(C2) 

(C3) 

(C4) 
(C5) 

Note that in the continuous limit where w\ — > w 2 so f\ 2 — > 1, we obtain the expected results that C 2 = C\ and 
D 2 = 0. It is interesting to note that it is the combination of the changes in w and c s that controls whether or 
not we have a non-trivial phase transition. 

We would like to calculate the power spectrum of the modes which are super-horizon at the end of inflation, ?/ e = 0, 
corresponding to x 3 (r) e ) <C 1. The transfer function is defined in Eq. (|37p . As described in the text we divide the 
modes of interest into three categories. The first category is defined by k > k\ 2 , the second category by ^23 < k < k\ 2 
and the third category by k < k 23 . 

First consider the first category, k > k\ 2 . For these modes we have £1(7712), £2(7712), £2(7723), £3(7723) S> 1. Using 
the large value approximation of the Hankel functions we obtain 



where we have defined 



fij 



sgn(l + 3w t ) (1 + Wj) c s 
sgn(l + 3wj) (1 + Wj) c s 



C 2 
C 3 



iDo 



fi 



1 



1 

iD 3 -e 



-fu + 1 

2i(x3(r]23)—Tr V3/2) 



" 2x 2 (ryi 2 )^- 1 / 2 



(/12 



^\gi(o:i(j)i2)-a!2(?7i2)-7r(i'i-i/2)/2) (C6) 

(C7) 



-C2 ^2(7723) 



1/2-1/2 



where 



1 
J 



2X3(^23)^- 1/2 



(1 + /aa) 
(1 - /as) 



; i(£c 2 ('j23)-a:3(';23)- , n'(^2-i'3)/2) 



1 - /23 1 - /: 



12 g2i(X2(»712)— Xa(fj2a)) 



1 + /23 1 + /l2 

1 + /23 1 ~ /l2 e 2i( £C2 (»7i2)-X2( t ;23)) 
l — /23 1 + /l2 



/23 
—I. 



1 



2i(a 2 ('7i2)-a;2(T723)) 



(C8) 



The last approximate equality is valid for the case in which w\ ~ w 3 ~ — 1 and as a result f\ 2 « 1 and /23 >> 1. 
Plugging these values of C 3 and D 3 in the formulas for the final curvature perturbation amplitude (|3"5j) and the 
transfer function (|3"T|) . we obtain Eqs. (|39l) and (|40p. 

Now consider the second category, /«23 < k < k% 2 . For these modes we have £1(7712), £2(^12) *C 1 while £2(7723), 
£3(7723) ~> 1. For this category, we obtain 



C 2 



£) 2e 2^ 2 _ _2^-^- 1 Cir(^ 1 )r(-^ 2 + i)e^ 2 



(C9) 



C 3 ~ i£>g^L e - 2< < B »- ,r, */ 2 > ~ ^C , 2£ 3 (7723) _ " 3+1/2 2;2(?723)" 2_1/2 e l(:C2( ' ,23) ~ a;3( ' ,23) "" (,y2 " ,y3)/2) , (CIO) 
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where 

H = (/23 + 1) + »(1 - f 2i ) e -^2(m^^2/2) _ _ ze -2 4 ( K2 (r, 23 )-3 7r i, 2 /2)^ _ (cu ) 

J" = (-/ 23 + l)+i(l + / 2 3)e- 2i(a;2(,)23) - 3 ^ 2/2) ~-X'. 

Again the second approximate equality holds for the case /23 S> 1. Now plugging these expressions of C 3 and D 3 into 
the formulas for the final curvature perturbation amplitude (1351) and the transfer function (f37|) . we obtain Eqs. (|42)) 
and (I43|) . 



2. Matter Dominated Intermediate Era 



Now consider the case where the intermediate non-inflationary stage is a matter-dominated universe, so w 2 — 
c s2 = 0- As described in the text for this case we work with the Bardeen potential $ while for the final curvature 
perturbation power spectrum we can switch to 1Z as given by Eq. (|53l) by using Eq. (|19[) . Our task is to find the 
expressions of A 3 and B 3 in Eq. (1531) in terms of A±, which is expressed in terms of G\ in Eq. (1 5 1) . 

Imposing the matching conditions (|18l) and (14"51) at rj — r\\ 2 and 77 = 7723 yields 



A 2 = - 



2wi Ax £1(7712 
5(1 +101) 



,1/1—1 



Xl(Vl2) ^(1) , , u „(1) / , u 

p iWi ^-2(^1(^12)) - ff^_i (a;i (7712)) 



B 2 = A ia :i( 77l2 ) l ' 1 - 1 ^ ) _i(xi(r ?12 ))-A 2! 



.4, 



8/ 



x 3 {r] 23 ) 



1-^3 



2^3(7723) \B 2 + A 2 



n 



+ fa I -2 W3 B 2 + (5 + 3w 3 )A 2 1^)1 H^>(x 3 (r, 23 )) 



B 3 = -^rX 3 (r] 23 ) 



1-^3 



2^3(7723) [b 2 +a 2 



H 2 3 
^12 
5\ 



r(2) 



(C12) 
(C13) 

(CM) 



2:: 



+ /? 3 -277J3B2 + (5 + 377J 3 )A 2 



"H 2 3 
"Hl2 



"Hl2 
5\ 



HSLl(X3(»&3)) 



(C15) 



As before, we evaluate A 3 and B 3 for the three different categories of the modes separately. For the first category, 
k > fci2, we have 



2aia : i(77i2)^- 1 / 2 [2i (xi ^ Ul/2) 
A *- B *- 5/3i(l+ W i) V7 6 

A 3 ~ -ie-^'^^^Bz^J^-BiX^) 



-V3+3/2 -i(x 3 (ri23)-nv3/2) 



For the second category, A23 < k < k\ 2l we have 



-i2" /1 77Ji 
3 + 5?«l ' 57r(l + 77)1 ) 



-B 2 



Air(i^-i), 



7T /2i 

4 V 7T 

Plugging these into the transfer function formula (|54|) yields Eqs. (|55|) and ([56 



A 3 ~ _ ie -2i(»3(%3)-^ 3/2)5 ^ 1 J ^ B 2 x 3 (r ]23 )-^ +3 / 2 e- l(x ^-^^ . 

4 V 7T 



(C16) 
(C17) 

(C18) 
(C19) 



Appendix D: Necessary relations for numeric calculations 



In order to have an efficient and accurate code, it is much better to change the variable from the conformal time to 
the number of e-folds, 

dn = Hdr). (Dl) 
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Using the Friedmann equation one has 



(l + 3w)dri) , 



H(ri) = H e exp 



where n e — n is the total number of e-folds from the time n until the end of inflation. 
The curvature perturbation obeys the equation of motion in terms of n, 



Similarly one has 



dn 2 



d 2 ® 
dn 2 



(» (1 _„ ) + £ lB((1+ „ )/ 4))<* + 



H 2 



-K = 0. 



- + 3(c 2 w ~w/2) 



d<I> 
dn 



2 1,2 



n 2 



3{w-c 2 w ) 



$ = 0. 



To model smooth transitions, we use the error function. For example the equation of state is given by 

w = w 1 + (-w-l + w 2 ) (1 + erf[di(n - rii)]) /2 + (-w 2 + w 3 ) (1 + er{[d 2 {n - (m + n 2 ))\) /2 , 
where d\ and d 2 determine the sharpness of the first and second phase transitions, respectively. 



(D2) 



(D3) 



(D4) 



(D5) 



References 



[1] E. Komatsu et al., "Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpreta- 
tion," larXiv:1001.4538l [astro-ph.CO]. 

[2] A. A. Starobinsky, "Spectrum of adiabatic perturbations in the universe when there are singularities in the inflation 
potential," JETP Lett. 55, 489 (1992) [Pisma Zh. Eksp. Teor. Fiz. 55, 477 (1992)]. 

[3] S. M. Leach, M. Sasaki, D. Wands and A. R. Liddle, "Enhancement of superhorizon scale inflationary curvature pertur- 
bations," Phys. Rev. D 64, 023512 (2001) astro-ph/0101406] . 

[4] J. A. Adams, B. Cresswell, R. Easther, "Inflationary perturbations from a potential with a step," Phys. Rev. D64, 123514 
(2001). [astro-ph/0102236] . 

[5] J. -O. Gong, "Breaking scale invariance from a singular inflaton potential," JCAP 0507, 015 (2005) astro-ph/0504383 . 
[6] M. Joy, V. Sahni, A. A. Starobinsky, "A New Universal Local Feature in the Inflationary Perturbation Spectrum," Phys. 

Rev. D77, 023514 (2008). [arXiv:071 1.15851 [astro-ph]] . 
[7] X. Chen, R. Ea sther and E. A. Lim, "Large Non-Gaussianities in Single Field Inflation," JCAP 0706, 023 (2007) 

|astro-ph/0611645] . 

[8] X. Chen, R. Easther and E. A. Lim, "Generation and Characterization of Large Non-Gaussianities in Single Field Inflation," 

JCAP 0804. 010 (2008) |arXiv:0801.3295l [astro-ph]]. 

[9] X. Chen, "Primordial Features as Evidence for Inflation," JCAP 1201, 038 (2012) |arXiv:1104.T323l [hep-th]]. 
[10] X. Chen, "Fingerprints of Primordial Universe Paradigms as Features in Density Perturbations," Phys. Lett. B 706, 111 

(2011) arXiv: 1106. 16351 [astro-ph.CO]] . 

[11] S. Pi and M. Sasaki, "Curvature Perturbation Spectrum in Two-field Inflation with a Turning Trajectory," arXiv: 1205.0161 
[hep-th]. 

[12] S. Hotchkiss and S. Sarkar, "Non-Gaussianity from violation of slow-roll in multiple inflation," JCAP 1005, 024 (2010) 

|arXiv:0910.3373l [astro-ph.CO]]. 
[13] F. Arroja, A. E. Romano and M. Sasaki, "Large and strong scale dependent bispectrum in single field inflation from a 

sharp feature in the mass," Phys. Rev. D 84, 123503 (2011) l arXiv: 1106.53841 [astro-ph.CO]]. 
[14] P. Adshead, C. Dvorkin, W. Hu and E. A. Lim, "Non-Gaussianity from Step Features in the Inflationary Potential," Phys. 

Rev. D 85, 023531 (2012) |arXiv:1110.3050l [astro-ph.CO]]. 
[15] A. A. Abolhasani, H. Firouzjahi, S. Khosravi and M. Sasaki, "Local Features with Large Spiky non-Gaussianities during 

Inflation." larXiv: 1204.37221 [astro-ph.CO], 

[16] F. Arroja and M. Sasaki, "Strong scale dependent bispectrum in the Starobinsky model of inflation," arXiv: 1204.6489 

[astro-ph.CO]. 

[17] A. E. Romano and M. Sasaki, "Effects of particle production during inflation," Phys. Rev. D 78, 103522 (2008) 
|arXiv:0809.5T42 [gr-qc]]. 

[18] D. Battefeld, T. Battefeld, H. Firouzjahi and N. Khosravi, "Brane Annihilations during Inflation," JCAP 1007, 009 (2010) 
|arXiv: 1004.14171 [hep-th]] . 



19 



[19] H. Firouzjahi and S. Khoeini-Moghaddam, "Fields Annihilation and Particles Creation in DBI inflation," JCAP 1102, 

012 (2011) [arXiv: 101 1.45 00 [hep-th]]. 
[20] D. Battefel d, T. Battefeld, J . T. Giblin, Jr. and E. K. Pease, "Observable Signatures of Inflaton Decays," JCAP 1102, 

024 (2011) [arXiv:1012. 13721 [astro-ph.CO]] . 
[21] N. Barnaby and Z. Huang, "Particle Production During Inflation: Observational Constraints and Signatures," Phys. Rev. 

D 80, 126018 (2009) |arXiv:0909.075T1 [astro-ph.CO]]. 
[22] N. Barnaby, "On Features and Nongaussianity from Inflationary Particle Production," Phys. Rev. D 82, 106009 (2010) 
arXiv: 1006.4615 [astro-ph.CO]]; N. Barnaby, "Nongaussianity from Particle Production During Inflation," Adv. Astron. 

2010, 156180 (2010) |arXiv: 1010.55071 [astro-ph.CO]]. 
[23] T. Biswas, A. Mazumdar and A. Shafieloo, "Wiggles in the cosmic microwave background radiation: echoes from non- 
singular cyclic-inflation," Phys. Rev. D 82, 123517 (2010) [arXiv:1003.3206l [hep-th]] . 
[24] E. Silverstein and A. Westphal, "Monodromy in the CMB: Gravity Waves and String Inflation," Phys. Rev. D 78, 106003 

(2008) [arXiv:0803.3085 l [hep-th]]. 
[25] R. Flauger, L. McAllister, E. Pajer, A. Westphal and G. Xu, "Oscillations in the CMB from Axion Monodromy Inflation," 

JCAP 1006, 009 (2010) |arXiv:0907.2916l [hep-th]] . 

[26] R. Flauger and E. Pajer, "Resonant Non-Gaussianity," JCAP 1101, 017 (2011) |arXiv:1002.0833l [hep-th]] . 

[27] R. Bean, X. Ch en, G. Hailu, S. -H. H. Tye and J. Xu, "Duality Cascade in Brane Inflation," JCAP 0803, 026 (2008) 

|arXiv:0802.049T1 [hep-th] ] . 

[28] M. Nakashima, R. Saito, Y. -i. Takamizu and J. 'i. Yokoyama, "The effect of varying sound velocity on primordial curvature 
perturbations," Prog. Theor. Phys. 125, 1035 (2011) |arXiv:1009.4394l [ astro-ph.CO]]. 

[29] V. Miranda, W. Hu and P. Adshead, "Warp Features in DBI Inflation," larXiv:1207.2T86l [astro-ph.CO]. 

[30] J. Emery, G. Tasinato and D. Wands, "Local non-Gaussianity from rapidly varying sound speeds," arXiv: 1203.6625 [hep- 
th]. 

[31] A. A. Abolhasani, H. Firouzjahi and M. Noorbala, "Effects of Variable Newton Constant During Inflation," arXiv: 1206.0903 
[astro-ph.CO]. 

[32] R. Allahverdi, A. Mazumdar and T. Multamaki, "Large tensor-to-scalar ratio and low scale inflation," arXiv:0712.2031 
[astro-ph] . 

[33] B. A. Bassett, S. Tsujikawa and D. Wands, "Inflation dynamics and reheating," Rev. Mod. Phys. 78, 537 (2006) 
|astro-ph/0507632| . 

[34] J. Garriga and V. F. Mukhanov, "Perturbations in k-inflation," Phys. Lett. B 458, 219 (1999) [hep-th/9904176], 

[35] I. Zaballa and M. Sasaki, "Boosted perturbations at the end of inflation," JCAP 1003, 002 (2010) |arXiv:0911.2069l 
[astro-ph.CO]]. 

[36] B. J. Carr, K. Kohri, Y. Sendouda and J. 'i. Yokoyama, "New cosmological constraints on primordial black holes," Phys. 

Rev. D 81, 104019 (2010) |arXiv:0912.5"297l [astro-ph.CO]]. 
[37] M. Drees and E. Erfani, "Running-Mass Inflation Model and Primordial Black Holes," JCAP 1104, 005 (2011) 

|arXiv: 1102.23401 [hep-ph]] . 

[38] M. Drees and E. Erfani, "Running Spectral Index and Formation of Primordial Black Hole in Single Field Inflation Models," 

JCAP 1201, 035 (2012) |arXiv: 1110.60521 [astro-ph.CO]]. 
[39] A. S. Josan, A. M. Green and K. A. Malik, "Generalised constraints on the curvature perturbation from primordial black 

holes," Phys. Rev. D 79, 103520 (2009) |arXiv:0903. 31841 [astr o-ph.CO]]. 

[40] A. D. Linde, "Hybrid inflation," Phys. Rev. D 49, 748 (1994) [arXiv:astro-ph /9307002] . 

[41] A. R. Liddle, D. H. Lyth, E. D. Stewart and D. Wands, "False vacuum inflation with Einstein gravity," Phys. Rev. D 49, 
6410 (1994) |arXiv:astro-ph/940101l] . 



